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ABSTRACT 

Polarization measurements from relativistic outflows are a valuable tool to probe the 
geometry of the emission region and the microphysics of the particle distribution. 
Indeed, the polarization level depends on: (i) the local magnetic field orientation, 
(ii) the geometry of the emitting region with respect to the line of sight, and (iii) the 
electron pitch-angle distribution. Here we consider optically thin synchrotron emission 
and we extend the theory of circular polarization from a point source to an extended 
radially expanding relativistic jet. We present numerical estimates for both linear 
and circular polarization in such systems. We consider different configurations of the 
magnetic field, spherical and jetted outflows, isotropic and anisotropic pitch-angle 
distributions, and outline the difficulty in obtaining the reported high level of circular 
polarization observed in the afterglow of GRB 121024A. We conclude that the origin of 
the observed polarization cannot be intrinsic to an optically thin synchrotron process, 
even when the electron pitch-angle distribution is extremely anisotropic. 

Key words: radiation mechanisms: non-thermal 


1 INTRODUCTION 

Circular polarization at the level of 0.6% has been 

recently detected in th e optical afterglow of GRB 121024A 
llWiersema et al.ll20l4 l. This is the first time that circular 
polarization is detected in GRB afterglow radiation. The 
same burst also shows linear polarization at a level of ~ 4%, 
implying a ratio p'^’-rc^phn ^ o.l5. Theoretical estimates of 
circular polarization in synchrotron radiation predict that 
both the absolute value of and its value relative to the 
amount of linear polarization p‘^"c^piin < 1/je, where 
ye is the random Lorentz factor of the radi ating electrons 
IILegg fc Westfoldlll968l : [Sazonov|[i969l . Il97^ 1. These calcu¬ 
lations have been performed in the frame of the fluid, and 
for a point-like region where the magnetic field has a given 
orientation. Based on these order of magnitude estimates 
and assuming that synchrotron radiation is dominating the 
optical afterglow of GRB 121024A, the measured values of 
both p=»’''= and jp^'-n large as compared to expec¬ 

tations, and their theoretical interpretation is challenging. 
For this burst, indeed, 7e is estimated to be of the order of 
10^ at the time when observations are performed, implying 
that observations are well in excess of the predicted value. 

A certain degree of linear polarization in GRB afterglow 
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radiation is instead expected, and has been ind eed measured 
in th e optical afterglow of several GRBs (see iGovino et al.l 
l2004l for a review). In all cases, the reported values of P**" 
is at the level of a few percent, with the noticeable excep¬ 
tion of GRB 120308A, where a higher level of polarization 
^pitri ^ 28%) has been measured in t he early (four minute s 
after the prompt) optical afterglow llMundell et al.l [201311 . 
Besides being a confirmation of the synchrotron nature of 
the afterglow emission, the level of polarization has been 
used to infer the properties of the magnetic field configura¬ 
tion and the geometry of the emission. The level of polariza¬ 
tion, indeed, depends on local conditions (as the orientation 
of the magnetic field and the properties of the electron popu¬ 
lation responsible for the emission) and on global conditions, 
as the geometry of the outflow. For this reason, polarization 
measurements are in general considered a valuable tool to 
learn about the physics of the source. 

If the magnetic field is fully tangled and has the same 
strength in each directions then the radiation is unpolarised, 
no matter the geometry of the outflow. A symmetry break 
in the configuration of the magnetic field must be intro¬ 
duced in order to explain the observed polarization. In the 
case of a spherical outflow, a completely tangled magnetic 
field in the plane of the shock (i.e. with no component in 
the direction perpendicular to the plane) would produce no 
net polarization, and some degree of coherence of magnetic 
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field lines must be invoked in o rder to break the symme¬ 
try (iGruzinov fc WajcmanI [T999l l . The situation is different 
in case of collimated outflows. In this case in fact, even a 
completely random magnetic field in the plane of the shock 
can give rise t o a net polarization, prov i ded that t he jet 
is seen off-axis (iGhisellirii fc Lazzatilll999l : ISar 1 119991) . This 
model predicts a change by 90° in the polarization position 
angle around the jet break time and a typical evolution of 
the degree of linear polarization. Different predictions on the 
evolution of p**" and its position angle are derived if an or¬ 
dered component and/or a compon ent perpendicular to the 
plane of the shock is considere d, llGranot fc Konigll l2003l : 
iGranotl l2003l : iNakar et al.l l2003l ) or if the magnetic field is 
significantly different in the parallel and perpendicula r di¬ 
rectio ns and sideways spreading of the jet is considered dSaril 
I 1999 I) . The temporal evolution of the linear polarization has 
also been proposed as a tool to reveal the structure of the 
jet, based on the fact that polarization curves for structured 
j ets are very diffe rent from those derived for uniform jets 
dRossi et al.ll2003 ). 

While in general observations of linear polarization can 
be explained within the current theoretical models, the 
detection of circular polarization was an unexpected re- 
sult. In the original pa per where the detection is presented 
dWiersema et al.ll2014l ). the authors show that it is very un¬ 
likely that its origin could be ascribed to plasma propagation 
effects or dust scattering, and argue in favour of an intrinsic 
origin. In particular, they invoke an anisotropic distribu¬ 
tion of the electron pitch angles cis a possible explanation. 
However, estimates are performed in the local frame of the 
fluid and do not take into account that the contribution of 
different emitting regions to the total polarization must be 
averaged, in order to derive the total (integrated over the 
emitting region) polarization in the frame of the observer. 
The relativistic motion of the outflow (that requires appli¬ 
cation of Lorentz transformations) and the integration over 
the (unresolved) image both play an important role in the 
determination of the final polarization detected by the ob¬ 
server and should be properly treated. 


In this work we derive the circular polariza¬ 
tion of synchrotron radiation from a relativistic out¬ 
flow. The amount of linear polarization for differ¬ 
ent magnetic field configurations and outflow ge ome- 
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Nakar et al.ll2003l: 

Rossi et al.ll2004|). However, since we are 


also interested in estimating the ratio j^ also 
report the equations to derive the linear polarization. 


The paper is organised as follows. In ^we provide an 
order of magnitude estimate and a schematic description of 
the origin of polarization in synchrotron radiation. Section 
^introduces the Stokes parameters in the comoving frame 
(i.e. the frame at rest with the source). In m the geomet¬ 
rical set-up we are considering is presented. Lorentz trans¬ 
formation to move from the comoving to the observer frame 
are presented in ^S] The equations to derive the total po¬ 
larization from an unresolved source are derived in ^ In 
El we estimate the circular and linear polarization for dif¬ 
ferent configurations. We discuss our results in comparison 
with polarization measurements of the optical afterglow of 




Figure 1. Upper panel; polarization of the radiation emitted by 
an electron with pitch angle a as detected by different observes. 
The angle represents the angle between the electron velocity /3e 
and the direction of the photon. Lower panel; polarization from 
a population of electrons with different pitch angles, as detected 
by one single observer, ip is the angle between the magnetic field 
and the direction of the observer. 


GRB121024A in section ^8] Finally, in 39]we summarise our 
conclusions. 


2 THE ORIGIN OF POLARIZATION IN 
SYNCHROTRON RADIATION 

In this section we discuss what is the origin of polarization 
in synchrotron radiation and why the level of circular polar¬ 
ization is expected to be small as compared to the level of 
linear polarization. This discussion is also aimed at under¬ 
standing which are the physical parameters the determine 
the expected level of polarization in the comoving frame of 
the fluid. 

We first consider emission from one single electron with 
a pitch angle a, gyrating in a uniform magnetic field B 
(Figure [T1 upper panel). Let tp be the angle between the 
electron velocity /3e and the direction of the photon n, 
and let consider three different observers detecting radiation 
from the same electron. In general, the observed radiation 
is elliptically polarised, with the axes of the ellipse parallel 
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and perpendicular to the projection of B on to the plane 
transverse to the direc tion of the photon n llWestfoldll 19591 : 
iLegg fc Westfoldlll968l ). The major axis of the polarization 
ellipse is perpendicular to the projection of B for |i/)| close 
to zero, but as |'!/)| increases, the form of the ellipse varies 
to a circle and then to an ellipse with major axis parallel to 
the projection of B. The direction is right-handed (RH) or 
left-handed (LH) according with ip ^ 0. The polarization is 
linear if -i/i = 0. We note a negligible amount of radiation is 
radiated at angles \ip\ > l/ 7 e = \/1 — /3|. 

Let us now consider a population of electrons with dif¬ 
ferent pitch angles, and a single observer, located at an an¬ 
gle ip from the magnetic field (Figure [TJ lower panel). As 
follows from the previous discussion, photons corresponding 
to Ip = 0 (i.e., emitted in the direction of the observer by 
electrons with a = ip) are linearly polarized. Photons for 
which Ip ^ 0 instead will be elliptically polarised. The ellip¬ 
tical polarization is the same for photons at ip and —ip, but 
with opposite direction of the rotation. 

The resulting polarization detected by the observer is 
obtained after integration over the contributions from all 
photons that reach the observer. First we note that only 
electrons with pitch angle ip—l/"/e < a < ip+l/js give a non 
negligible contribution to the emission detected by the ob¬ 
server. These pitch angles correspond to photons with very 
small Ip. As noticed before, small ip correspond to polariza¬ 
tion ellipses with major axes perpendicular to the projec¬ 
tion of B on the plane perpendicular to the observer. This 
implies that when the contribution from all these ellipses 
with all major axes perpendicular to the projection of B is 
summed up the final level of integrated lineal polarization is 
still large (around 70%). On the contrary, the integrated cir¬ 
cular polarization, is largely reduced by integration. Indeed, 
the contribution from photons with ip > 0 is partially can¬ 
celled by the contribution from photons with ip < 0, since 
different signs of ip correspond to different directions of the 
rotation. This cancellation suppresses the final amount of 
circular polarization, which is expected to be small as com¬ 
pared to the linear polarization. The level of cancellation 
depends on how many photons have positive and negative 
Ip-, if the number of photons with angle ip is equal to the 
number of photons with angle —ip no net circular polariza¬ 
tion is detected. The number of photons at different angles 
depends on the number of electrons at different pith angles. 
This is the reason why, besides the angle ip, the properties 
of the polarization obtained after averaging over the con¬ 
tributions from different electrons depend also on the pitch 
angle distribution and on the shape of the electron energy 
spectrum. 

In the following we consider an electron distribution 
that has a power-law dependence on the energy, the num¬ 
ber of electrons with Lorentz factor 7 ^ per unit solid angle 
around the pitch angle a is: 

N{a,'f^) cxY{a)’y~\ ( 1 ) 

With the index <5 we refer to the electron distribution after 
radiative cooling has modified the initial injected distribu¬ 
tion. Hence this distribution can be different than the in¬ 
jected electron spectrum Ninjpe) oc . For an uncooled 
distribution S = p, while for cooled electrons <5 = p + 1 or 
5 = 2, depending on the value of 7 e as compared to the 
minimum and the cooling Lorentz factors. 


The function Y{a) is defined as: 


Yia) 


dN{E, a) 47 r 

dn NeiE)' 


( 2 ) 


where Ns{E) is the total number density of electrons with 
energy E. Following this definition, Y = 1 for an isotropic 
pitch angle distribution. 

The influence of the pitch angle distribution is clear 
(see Figure [!}: if the number of electrons with pitch angle 
Q = ip+ip is different from the number of electrons with pitch 
angle a = ip — ip, the contributions to the total polarization 
from electrons at ip and —ip do not cancel out. For example, 
for the case of isotropic pitch angle distribution, if (p < 90° 
the number of electrons with -i/i > 0 is larger than the number 
of electrons with ip < 0 and the total polarization is negative 
(as in the case in Figure [T]). Viceversa, if p > 90° the total 
polarization is positive. 

The role played by the power-law index <5 of the Lorentz 
factor distribution can be understood by first noticing that 
for a given fixed frequency 7 ^ oc 1/ sin a. Electrons with dif¬ 
ferent pitch angles (as the ones represented in Figure[T]) must 
then have different Lorentz factors in order to radiate pho¬ 
tons at the relevant frequency i/. For this reason the overall 
polarization is also determined by the number of electrons 
at different 7 e, i.e. from <5. 


3 LOCAL STOKES PARAMETERS 

The Stokes parameters I, Q, U, V describe the polarization 
state of the electromagnetic radiation. They represent ac¬ 
tual intensities, and more specifically, linear combinations of 
intensities measured in orthogonal polarizations directions. 
In terms of the parameters of the polarization ellipse (that 
has been qualitatively described in Section 2), the Stokes 
parameters for a monochromatic wave are: 

I={a^ + P), (3) 

Q = (a^ — P) cos2A, 

U = [P — P) sin 2A, 

V = ±2a6 

where a and b refer to the semi-major and semi-minor axes of 
the ellipse, and A is the tilt angle of the ellipse with respect 
to some reference direction. The parameter I describes the 
total intensity of the beam, Q describes the excess of linearly 
horizontally polarized light over linearly vertically polarized 
light, U describes the excess of linear -1-45° polarized light 
over linear —45° polarized light, and the fourth parameter 
V describes the excess of right-circularly polarized light over 
left-circularly polarized light. 

As explained before, the major and minor axis of the po¬ 
larization ellipse (and then, the Stokes parameters) depend 
on the angle ip between the photon direction and the elec¬ 
tron velocity. Then integration over all the radiation must 
be performed in order to find the Stokes parameters coming 
from a population of electrons, located at some positions of 
the source. In the next sections we give the equations for the 
Stokes parameters averaged over the electron population, in 
the comoving frame of the source. We add the subscript ’0’ 
in order to distinguish these local Stokes parameters (i.e. 
estimated for a point-like region) from those averaged over 
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the spatially extended emission region. Since the linear and 
circular polarization are a combination of the normalised 
Stokes parameters (i.e. relative to the total intensity of the 
beam, described by the first parameter Jo), in the following 
we give te equations for Uq/Iq, Qo/Iq, and Vq/Iq. 


3.1 Linear polarization: the Stokes parameters 

Qo/Io and Uo/Io 

The linear polarization P^aa: from an isotropic electron dis- 
tribution, and for a given d irection of the magnetic field is 
dRvbicki fc Lightman [l986l l : 


Pr. = (<5+1)/(^ + 7/3). 


(4) 


This quantity can also be expressed in terms of the spectral 
index s of the radiation spectrum (defined by dF^/dE oc 
Recalling that the spectral index of radiation emitted 
by a population of electrons with a power-law energy spec¬ 
trum {dNe/dE oc E~^) is given bys = {S — l)/2, we obtain 
Pmax = (s -I- l)/(s -I- 5/3). 

For configurations where locally the magnetic field is or¬ 
dered (like uniform in the plane of the shock, or perpendic¬ 
ular to the plane of the shock), the local linear polarization 
is: 

Po‘" = PmL (orderd), (5) 


and does not depend on the position. For typical values of 
S, the linear polarization estimated from a point-like region 
is very high (as predicted in the quantitative description in 
Section [2]), and is around ~ 60-70%. 

For configurations where the magnetic field is locally 
tangled, the polarization is obtained after averaging over all 
possible directions of the magnetic field. In the case of a 
perfectly random magnetic field i n the plane of the shock 
l|Sarilll99fll : lOranot fc Konig]ll2003l l : 

phn ^ (random in the plane). (6) 

1 + cos^ 8 ' 

The angle 9' is the angle between the direction of the photon 
(in the comoving frame) and the normal to the plane where 
the magnetic field lies. In this case the local polarization 
depends on the position of the observer as compared to the 
plane of the shock. 

The linear polarization is described by Qo/Io and Uo/Io- 
Qo/Io = Po"cos{2ep) (7) 

Uo/Io = Po” sm{2ep) 

The angle 9p is the position angle and gives the orienta¬ 
tion of the linear polarization in some coordinate system. In 
configurations with ordered magnetic field, even if the lo¬ 
cal polarization is independent on the location of the source 
(see equation[5]), the polarization angle is not, and when the 
contribution from different regions is summed partial cancel¬ 
lation can arise, depending on the geometry of the system. 


3.2 Circular polarization: the Stokes parameter 

Vo/Io 

The Stokes parameter Vo/Io describes the level of circular 
polarization. Its deri vation for a power-la w distribution of 
electrons is given by (ISazonovlll969l . 0*97211 : 



Figure 2. Representation of the geometrical set-up. In the frame 
at rest with the fluid, the photon is emitted at an angle 9' from 
the direction of the fluid velocity /?> but in the laboratory frame it 
is beamed in the direction of the observer. The emitted radiation 
is then mapped onto the plane of the sky, where we introduce 
a system of polar coordinates (r, iji). The coordinate r is equal 
to unity for a region located a.t 9 = l/F. The contribution from 
radiation at r > 1 is negligible. 


j^circ _ 

^0 = 


lo 


A^{5 + i){5 + 2)r{^)r{^) 

35(5 + 7/3) r(2^)r(2^) 


X 


(^cot ip + 


^_1 dY{p) \ 

5 + 2 Y{ip) dp ) 


1 

7e' 


( 8 ) 


This equation is valid provided that | '^^dUp'^ I < 7e . The 

circular polarization can also be expressed in terms of the 
linear polarization Pq™ = Pj/ax- 


j^circ 

-^0 


(^cot p + 


^_1 dY{p) \ 

S + 2Y(p) dp }' 


(9) 


where we have neglected factors order unity. 

In the standard case of isotropic pitch angle distribution 
(where the second term of the sum in brackets vanishes) the 
ratio between circular and linear polarization is of the or¬ 
der of l/ 7 e -C 1, and since P/////xx ~ 0.7 the value of 
is also of the order l/ye -C 1. Locally the circular polariza¬ 
tion is strongly suppressed and is of the order l/7e -C 1. As 
anticipated in Section [2l the reason for this suppression is 
the cancellation between electrons contributing with positive 
and negative polarization. A high level of circular polariza¬ 
tion can be obtained only from an electron population with 
a high level of pitch angle anisotropy. However, this is not 
the end of the story. Equation in fact cannot be compared 
directly to the observations. Comparison with observations 
requires that we first transform from the comoving to the 
observer frame and average polarization coming from differ¬ 
ent regions, as it will become clear in the next sections. 


4 THE GEOMETRICAL SET-UP 

In the previous section we have discussed polarization in the 
comoving frame, coming from electrons located in a specific 
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point-like region (i.e. a region where the magnetic field has 
a given orientation). However, we are interested in deriv¬ 
ing polarization from a blast wave expanding at relativistic 
velocity H 3> 1. This has two implications: i) we have a 
spatially extended emitting region that, as we will discuss, 
is characterised in each point by different local conditions, 
and ii) the system discussed in the previous section is now 
moving at relativistic speed with respect to the laboratory 
frame. In what follows, we describe the geometrical set-up 
that we are considering. A sketch is depicted in Figure [2] 
In this sketch the blast wave is represented as a sphere (or, 
equivalently a jet with jet opening angle Oj^t l/F and 
edges that are not visible to the observer). However, the 
procedure that we will present is also valid for jets with 
9jst < l/F and/or jets seen off axis. 

In the comoving frame, the direction of the photon is 
identified by the angle 8', measured from the direction of 
the bulk motion f3 = — l/F. Because of beaming, in 

the observer frame this photon is boosted in the direction of 
the observer, and makes an angle 6 with the fluid velocity (3 
(see Figure [2]). The relation between the two angles in the 
two different frames is given by a Lorentz transformation, 
and is described in Section [5l Besides the direction of the 
photon in the comoving frame, we also need the direction of 
the magnetic field in order to compute the polarization. This 
depends on the configuration of the field. However, it can be 
immediately realised that both the field component in the 
plane of the shock (i.e., tangential to the sphere surface) 
and the field component perpendicular to the plane of the 
shock (i.e., parallel to the radial velocity /3) have different 
orientations at each point of the emitting surface. From this 
discussion it is clear that the observer will detect different 
polarization from different regions. 

The procedure to compute the resulting polarization de¬ 
tected by the observer (i.e., the quantity that must be com¬ 
pared to the observational measurements), is the following. 
First, one needs to estimate the local (i.e. in a small region) 
polarization (that we have described in terms of local Stokes 
parameters) in the comoving frame of the source (Section[^. 
Each different region will be characterised by a different po¬ 
larization status. Then one needs to apply Lorentz transfor¬ 
mation to move from the comoving frame to the observer 
frame (Section [^l. Since the region from which the radia¬ 
tion originates is spatially extended, these different regions 
and their local polarization status are mapped into a plane 
perpendicular to the line of sight, i.e. the plane of the sky, 
and result in a polarization map (see Figure [2]). The emis¬ 
sion from GRB afterglows, however. Is spatially unresolved, 
so that the detected polarization comes from the integration 
over this polarization map. To estimate this integrated po¬ 
larization we first need to average each Stokes parameter. 
The average over the plane of the sky is, of course, weighted 
by the intensity of the emission. Finally, from the averaged 
Stokes parameters, the total integrated polarization is ob¬ 
tained (Section inj. 

In order to describe the local Stokes parameters in each 
point of the plane of the sky, we introduce a system of polar 
coordinates (r, <f)) on this plane. The center of this coordi¬ 
nate system is located at the intersection between the plane 
of the sky and the direction connecting the observer and the 
central engine. We define r = sin OF ~ OF (for r F), 
so that r = 1 corresponds to radiation coming from a re¬ 


gion located at sinS = 1/F. This is the visible region for 
the observer: electrons travelling at angles 8 > 1/F do not 
significantly contribute to the radiation received by this ob¬ 
server. We will also use the coordinate y in place of r, defined 
as y = (dF)^ = P. We call polarization map the mapping of 
the polarization (in the lab frame) coming from different lo¬ 
cations of the surface. We will give examples of polarization 
maps in section [T] 


5 THE LORENTZ TRANSFORMATIONS 

In order to explicitly estimate the local Stokes parameters, 
we first need to express the angles 8p and ip and the Lorentz 
factor 7 e as a function of quantities in the observer frame. 
We apply these transformations to each point of the emitting 
region and give the relevant parameters 8p, ip, and 7e as a 
function of the polar coordinate system (r, />) introduced in 
the plane of the sky. 

We are interested in those photons that, after Lorentz 
boosting, travel towards the observer. These photons are 
those emitted at an angle 9 from the direction of the fluid 
velocity j3. In the local frame (primed quantities), this angle 
Is: 

rii-"pcos8Y 

where T> is the Doppler factor: 


F(1 — /3 cos9) 1 + y 

This equation (where we have used the approximation 
cos 6 ~ 1 — 9^ /2) shows that the Doppler factor depends 
on the position on the map. This is due to the fact that we 
are considering a spherical emitting surface (or a part of a 
sphere, in case of jetted outflows), where the fluid moves in 
radial directions. This implies that different parts of the fluid 
move at different angles 9 with respect to the line of sight. 
Since sin 0 ~ 0 we can write the angle 9' (that identifies the 
direction of the photon in the local frame) as a function of 
the coordinates in the laboratory frame: 

sln6i' = -^^. (12) 

Now to estimate 9p we proceed as follows. First we need 
to identify in each point the direction of the magnetic field. 
This will depend on the configuration considered. Once the 
configuration has been chosen, the polarization vector in 
the local frame can be derived from the vectorial product 
between the magnetic field vector and the direction of the 
photon. Now that the polarization vector has been identi¬ 
fied, it is necessary to move back to the laboratory frame and 
describe the orientation of the polarization vector in our co¬ 
ordinate system. We give the polarization angle for the three 
different magnetic field configurations that we explore in the 
following: i) uniform in the plane of the shock, li) random in 
the plane of the shock, and iii) radial, i.e., perpendicular to 
the plane of the shock. For case i) we consider a magnetic 
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field parallel to the axis (/> = 0: 


6 THE TOTAL POLARIZATION 


Op 


+ arctan ( - - - cot 6] 

\^ + y J 

Op = (jj 

Op = (j) + 90° 


(uniform), 

(13) 

(random), 

(14) 

(radial). 

(15) 


These equations for Op, together with equations [7] and |4] al¬ 
lows us to compute the normalised Stokes parameters Qo//o 
and Uq/Iq in the laboratory frame at each position on the 
plane of the sky. Note that is given by equation [5] for lo¬ 
cally ordered magnetic fields (namely, uniform in the plane 
of the shock or radial), while is given by equation |S] for a 
random magnetic field. In this last case, using equation 1121 
we can write the local linear polarization as a function of 
the coordinate y: Pq"" = x i//(l + y^). 

To derive Vo/Io as a function of the position on the 
plane of the sky, we need to calculate and 'ye{y,<j)). 

To derive 7 e we first note that the observed frequency is 
1 /°*’® = i//(l + z) = v'D/il z), where the Doppler factor 
T> is to be applied to transform from the comoving frame 
frequency iz' to the progenitor frame frequency i/. 

The comoving frequency v' is given be synchrotron fre¬ 
quency formula: 


, 3uh sin a 2 
^ -7e, 


(16) 


where vh = eBI {2TTmsc). The pitch angle a can be taken 
equal to ip, since in deriving equation [S] the average on the 
electrons with slightly different pitch angles (in the cone 
l/7e) has been already performed. By inverting equation 1 161 
and writing everything in terms of the observed frequency 
we obtain the Lorentz factor of the electrons that radiate at 
the relevant frequency: 


2 u°^^l + z){l + y) 

Suh sin pP 


(17) 


Equation [T7] shows that the contribution to the emission at 
a given observed frequency (at which polarization mea¬ 
surements are performed) comes from electrons with differ¬ 
ent 7 e. This is due to the fact that in order to move from 
the observer to the comoving frame we need to account for 
the Doppler factor T>{y) and that the relevant pitch angle 
of those electrons which are radiating in the direction of the 
observer also varies. 

Finally we need to estimate p. This is the angle between 
the photon direction and the field lines, and it will depend on 
the configuration of the magnetic field. For a given direction 
of the magnetic field, the angle p is obtained by estimating 
the scalar product (in the local frame) between the photon 
direction and the magnetic field. We give its value for a 
uniform magnetic field in the plane of the shock and for a 
radial magnetic field (in a random magnetic field the local 
circular polarization is always zero): 


sm p = 



2 

cos^ 4> + sin^ (j} 


sin p = 


1 + y 


(uniform). 

(18) 

(radial). 

(19) 


To derive the total polarization we integrate each Stokes 
parameter separately over its map on the plane of the sky, 
after the contribution coming from different regions has been 
weighted by the intensity. Namely, the value of the Stokes 
parameter in each point of the map is weighted by the ratio 
between the contribution of that region to the total flux and 
the total flux. 

We define as the specific flux per unit surface 

dS = d(j>rdr = d(l>dy/2 coming from a region located at 
(y, (f)). Its value is proportional to: oc I^dQ oc ■ 

Recalling that the specific intensity /(,, in the comoving 
frame of the fluid is proportional to the number of electrons 
pointing in the direction of the observer we obtain: 

oc Y{p)V^{iz')-‘‘{Bsmpy+^ oc Y{p)V^+‘’{Bsinpy+^ 
cto 

oc T(yj)(l-f y)-(^+“)(sin(^)'’+\ (20) 


The total flux is the integral of equation (p0)l on the 
region of the map that corresponds to the emitting surface 
dS: 

F. = j I Yipyi+y)-^^+‘\smpy+^d<pdy/2. 

( 21 ) 

The proportionality constants in equations [20] and [21] are 
the same, and we don’t need to explicitly write them, since 
we only need the ratio between the partial flux at a given 
position and the total integrated flux Fv Integra¬ 

tion must be performed at least up to 1 / > 1, since this is the 
region where most of the emissivity come from. Integration 
up to la rger distances can introduce differences at the level 
of 10% llGranotll2003h . 

If the outflow instead of being spherical is collimated 
into a jet, in the plane of the sky the jet can be described 
in terms of position of its axis (r^et, 4>jet) and opening angle 
Ojst- Note that Vjst = O^i^wP, where 0-aie-w is the viewing 
angle, i.e. the angle between the line of sight and the jet 
axis. For on-axis jets, r^et = 0 and = 0, and the edges of 
the jet are described by = OjetP- Outside this region the 
flux vanishes. More generally, for off-axis jets Fi,{y, <()) = 0 in 
the region of the plane that satisfies the condition: (r cos4> — 
r j et cos (t>jet)^ -I- (rsin<() — rjet sin (ftjet)^ > OjetF. 

Taking as an example Uq/Iq (but the same equations 
can be applied to the other normalised Stokes parameters) 
the flux-weighted parameter is given by: 



TT dF^iyp) 

% ^ y(v2)(l + y)-<^+^nsiny>)^+^ 

lo f Y(p)(l + y)-(^+‘’i(sinp)‘’+^dS 


( 22 ) 


The subscript fw stands for ’flux-weighted’. We recall that 
the normalised Stokes parameters are also a function of 
the position (y, tfi). The integration of the quantity in equa¬ 
tion 1221 over all the surface will give the Stokes parameters 
for an unresolved source: 


U_ 

I 


(UA f^xi^d^dy 

J f-w f ^Fid^d(j)dy 


(23) 
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Polarization in ultra-relativistic synchrotron sources 7 


The total linear polarization is given by 



The total circular polarization is given by: 


V 

I 




(24) 


(25) 


7 RESULTS: SOME EXAMPLES 

Using the above formalism we discuss now polarization es¬ 
timates for several configurations. For the geometry of the 
outflow we consider three different cases: 

• a sphere (or equivalently a jet with jet opening angle 
0jet » l/F), 

• a narrow jet with 9jet = 1/F, 

• a narrow jet with 9jet = 1/3F. 

For the jetted geometries, the results depend on the viewing 
angle 9view between the jet axis and the line of sight. 

For the pitch angle distribution we consider: 

• the standard isotropic pitch angle distribution (corre¬ 
sponding to Y{a) = 1), 

• anisotropic distributions. 

We consider configurations of the magnetic field where 
the symmetry is broken, since completely tangled magnetic 
helds do not give rise to net polarization. We consider: 

• a random filed in the plane of the shock (i.e., with no 
component in the radial direction), 

• a radial magnetic field (i.e., with no component in the 
plane of the shock), 

• a uniform magnetic field in the plane of the shock (i.e., 
the component in one direction is much larger than the com¬ 
ponent in the orthogonal direction). 

We divide the presentation of the results into two 
sections: isotropic and anisotropic pitch angle distribu¬ 
tion. We will perform calculations in the spectral range 
fm < V < Vc, i.e., in the range of frequency between the 
injection frequency Um and the cooling frequency Ve- In 
this Ccse the slope of the election distribution is equal to 
the slope of the injected distribution, since electrons are 
in slow cooling regime: 5 = p. We chose the value p = 2.5, 
that implies a spectral index s = (p — l)/2 = 0.75 and 


T.l Isotropic pitch angle distribution 


In case of isotropic pitch angle distribution, estimates 
of linear polarization have been already performed for 
all the magnetic field geometries that we are consider- 


ing 

Gruzinov & Waxmanll 19991: iGhisellini & Lazzatilll999l: 


19991: Granot & Koniell 20031: Granotll2003l: Nakar et al.l 

2oo; 

: Rossi et al.| 2004). 


We report the results for completeness and for compar¬ 
ison with our results on circular polarization, since we are 
also interested in estimating the ratio 


If the pitch angle distribution is isotropic then dY/dtp = 
0, and the equation for simply reduces to: 

= (26) 
7e 

The function /(5) depends very weakly on 5, and assumes 
the value /(5) ~ 2 in the range & = [2.1 — 2.8]. p and 7 e 
as a function of the position on the map can be derived 
through equation 1181 191 and 1 171 However, the estimate of 7 e 
depend on P and B. In order to perform calculations which 
are independent on the specific values of these two param¬ 
eters, Instead of estimating p§"’^ we estimate the quantity 
Pq"^ X ')e,min, where ye,min is the minimum random Lorentz 
factor of those electrons contributing to the emission at fre¬ 
quency u. As it will become evident later, x ye,min 

does not depend on P and B. The relevant value of ye,min 
can be estimated from the equations that we will give and 
by deriving the values of P and B from standard afterglow 
equations. Note that, since the local circular polarization is 
of the order of l/7e, the total polarization will always be of 
the same order or smaller, no matter the geometry of the 
emission region (e.g., spherical or conical, top hat or struc¬ 
tured jet, with uniform or angle dependent luminosity). 

7.1.1 A uniform magnetic field in the plane of the shock 

We now consider a uniform magnetic field in the plane of 
the shock. The equations given below and the polarization 
maps are derived for a horizontally oriented (from left to 
right) magnetic field. In this configuration the angle p be¬ 
tween the field and the photon (estimated in the comov¬ 
ing frame) is given by equation (TS] The Lorentz factor ye 
of electrons radiating at frequency v reaches its minimum 
value when y = 0, that is for electrons that move towards 
the observer (i.e., 0 = 0). In this case p = 90° and sin p = 1 
(see equation [17]): 


I 7 [ship . . 

From this last equation and eq. [26]we obtain: 

X y.,min = -if^COtp f{S)P^:,. (28) 

V i + y 

Figure [3] shows the polarization maps for the case of a 
uniform field. The upper panel is the map of x ye,min. 
Negative and positive values refer to the rotation direction 
of the polarization. For tf = ±90° the polarization is zero, 
since those photons that reach the observer are always per¬ 
pendicular to the magnetic field {p = 90°, see Figure [T]). If 
f) = Q (or 180°) and j/ = 1 then p = Q, 180°, and equation 
[8] can not be applied. When —90° < < 90°, the angle be¬ 

tween the magnetic field and the photon is > 90° and the 
polarization is positive, while it is negative in the opposite 
case (see the discussion in Section [2|). 

The central panel shows the map of the flux-weighted 
circular polarization (Po'’^°X7e,min)/u; estimated from equa¬ 
tion 1221 where Uo/Io in this case must be replaced by 

ppirc y 
0 ie,min ■ 

Because of left-right symmetry, in the case of a sphere, 
or a jet seen on-axis, or an off-axis jet with axis located at 
(f> = ±90° the total integrated polarization = Q, i.e., in 
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Figure 3. Polarization levels for the case of an isotropic pitch an¬ 
gle distribution and a uniform magnetic field in the plane of the 
shock. The coordinate system {x,y) = (r cos 0 , r sin <;i) adopted 
to build the maps is shown in Fig. |2] Note that r = 1 corre¬ 
sponds to an angular size l/F. Upper panel; shown is the max¬ 
imum local circular polarization ^ •yg Middle panel: 

flux-weighted maximum local circular polarization per unit sur¬ 
face. Bottom panel; flux weighted local linear polarization per unit 
surface. The dotted circles show the position of the jet needed to 
maximise the total circular polarization for a jet with dj^t = 1/F 
and 0jef = l/3r. 


this magnetic field configuration the total circular polariza¬ 
tion vanishes regardless of what the local values of circular 
polarization are. The integrated polarization can differ from 
zero only for an off-a:xis jet, and its value depends on Ojet and 
on the position of the jet a:xis with respect to the observer. 
We study the case of OjetT = 1 and OjetP — 1/3. polar¬ 
ization maps are similar to those obtained for the spherical 
case, except for the fact that outside the region of the jet 
the flux is zero. We derived the total polarization for each 
different position of the jet axis in the plane of the sky and 
report in table [T] its minimum and maximum value and the 
position at which these values are reached. The minimum 
polarization is zero, as previously discussed. The maximum 
polarization instead is reached when 0 = 0,180° and r — 0.3 
(r = 1) for the narrower (larger) jet. The position of the jet 
needed to observe the maximum polarization is shown in 
Figure [3] for the two cases Oj^tP = 1 and OjatP = 1/3 with 
a dotted circle. 

The bottom panel of Figure [3] shows the map of the 
flux-weighted linear polarization (before weighting for the 
flux the linear polarization is the same in each point of the 
map and is equal to P^2x)- The total polarization in the 
spherical case is P*™ = 0.61. Values of P*™ at different 
frequencies, for different values of p and for integration over 
Vmax = 1 and ymax S> 1 are reported in iGranot ll2003l ') . 

7.1.2 A random magnetic field in the plane of the shock 

The local circular polarization always vanishes for an 
isotropic pitch angle distribution. The linear polarization 
locally can be very high (see equation O, but due to the 
symmetry in this configuration, the total polarization van¬ 
ishes for spherical blast waves and for wide jets. Net total 
linear polarization only arises for jets where the edges are 
visible. 

7.1.3 A radial magnetic field 

In a radial magnetic field (i.e., perpendicular to the plane of 
the shock) B is parallel to the direction of the fluid velocity 
/3, and ifi = 6', leading to equation 1191 

In this configuration the angle p does not depend on 
0. The minimum electron Lorentz factor is reached when 
y = 1/3. The circular polarization (multiplied by 'ymin) is: 

Po"" X 7e,min = “ COt Ifi f {S) ■ (29) 

The map of this geometry is shown in the upper panel of 
Figure [J] while the middle panel shows the same quan¬ 
tity weighted by the flux. For y < 1 (i.e., the region that 
mainly contributes to the observed flux) the angle ip is al¬ 
ways smaller than 90° and the polarization is negative. Its 
contribution to the total polarization is not canceled by the 
positive contribution coming from the region at y > 1, and 
the total polarization is = —0.70/7e,min. Jetted ge¬ 

ometries can help to increase a bit this value (a factor of 
2-3), but since locally the circular polarization is small, no 
large polarization can be reached by integration over the 
map (see table [T|). 

Locally, the linear polarization is equal to Pj^ax — 70%, 
but when integrated over all the emission region the total 
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Figure 4. Polarization levels for the case of an isotropic 
pitch angle distribution and a radial magnetic field perpendic¬ 
ular to the plane of the shock. The coordinate system {x,y) = 
(r cos cj>, r sin 0) adopted to build the maps is shown in Fig. [2] 
Note that r = 1 corresponds to an angular size l/F. Upper panel: 
shown is the maximum local circular polarization ^ . 

Middle panel: flux-weighted maximum local circular polarization 
per unit surface. Bottom panel: flux-weighted local linear polar¬ 
ization per unit surface. 



Figure 5. Schematic representation of the cancellation arising 
from different portions of the emitting region when an anisotropic 
pitch angle distribution is considered. In this example it is as¬ 
sumed that only electrons with pitch angles between 30° and 
60° (shaded region) exist. The solid arrows show the edges of 
the distribution, while the dashed arrows show the magnetic field 
(assumed to be uniform in the plane of the shock). From region 
T’ the observer detects positive circular polarization, whose value 
can be high, due to the lack of electrons that would contribute 
with negative polarization. Conversely, from region ’2’ only elec¬ 
trons that contribute with a negative polarization are present. 
Outer regions do not contribute to the flux, since they are not 
emitting photons in the direction of the observer. 


linear polarization is null: P**" = 0. This can be understood 
by looking at the orientation of the polarization vector in 
Figure^] (bottom panel). If jetted geometries are considered, 
then P**" can differ from zero and reach values as high as 
40% for very narrow jets. 

From these considerations we conclude that for a radial 
magnetic field the ratio from l/7e to in- 

hnitive. However, the absolute value of p°®’’° is still limited 
to be < l/7e. 


7.2 Anisotropic pitch angle distribution 

For an isotropic pitch angle distribution, a large circular po¬ 
larization (i.e., in excess of l/7e) cannot be achieved. This 
statement is true not only for the total integrated polariza¬ 
tion, but also for a point-like region. In fact, as shown in 
Figure [T] radiation originating from a point-like region is 
emitted by electrons with different pitch angles, and partial 
cancellation suppresses the polarization. However, if elec¬ 
trons are missing for some pitch angles (namely, if the an¬ 
gle distribution is anisotropic) then the polarization from a 
point-like region can be higher. For example (see Figure [T]), 
if all the electrons with pitch angle a < ip are missing, then 
the contribution from photons with ip > 0 is not suppressed 
by contribution from photons with ip < 0. Then the con¬ 
dition for high Pp'’'° is that the observer is located just at 
the edge of the pitch angle distribution and that the level of 
anisotropy is very large. However, this is true if we consider 
the emission coming only from a specific location. Cancella¬ 
tion can still arise due to contribution from radiation emit¬ 
ted at different locations. 

Consider two different point-like regions as in Figure [5] 
and imagine that the magnetic field is uniform in the plane 
of the shock (dashed arrows), and that electrons can only 
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UNIFORM B 

RANDOM B 

RADIAL B 

piin 

pcirc y ^ 

-* ie,rmn 

piin 

pun 

pcirc y ^ 

-* ^ ie,rmn 

Sphere 

0.61 

0 

0 

0 

0.70 

II 

min 

0.61 

[4) = 0,180°;r = 1) 

0 

(</ = ±90°) 

II o 

o 

0 

(r = 0) 

0.81 
(r = 1) 

max 

0.66 

(0.1 <r < 0.85) 

0.66 

{(p = 0,180°; r = 1) 

0.10 

0.20 
(r = 1) 

1.0 

(0.49 <r < 0.63) 

Qj - i/sr 

min 

0.71 

0 

{<p = ±90°) 

II o 

o 

0 

(r = 0) 

1.6 

(r = 0.33) 

max 

0.72 

0.64 

((/> = 0,180°; r = 0.33) 

0.10 

0.40 

(r = 0.33) 

2.4 

(r = 0) 


Table 1. Isotropic pitch angle distribution; total linear and circular polarization for different geometries of the emitting region and 
different configurations of the magnetic field. The reported values for the circular polarization are multiplied by 7e.min- For jetted 
geometers the polarization depends on the position (r, ip) of the jet axis with respect to the observer, located at r = 0. We list the 
maximum and minimum polarization that can be reached, and report the corresponding position of the jet. When r and/or (p are not 
specified, it means that the value is independent on that parameter. In a random magnetic field circular polarization is always zero. 


have pitch angles in the range 30° — 60° (shaded areas). 
With reference to the point-like region T’, the observer is 
located at one of the edges of the distribution (solid arrows), 
so that the net polarization is high (due to the high deriva¬ 
tive of the distribution) and is positive. However, when we 
consider region ’2’, the observer now sees the opposite edge 
of the distribution and the net polarization is still high, but 
negative. If we consider intermediate regions between ’1’ and 
’2’ the polarization varies from large positive values down 
to zero, and then becomes negative with increasing abso¬ 
lute values, and reaches the maximum negative value at the 
location denominated ’2’. Outside the region ’1’ to ’2’ the 
amount of radiation that reaches the observer is negligible. 
The only way to avoid cancellation from the positive and 
negative regions is to consider a jetted geometry. If the jet 
is large enough that all the region from T’ to ’2’ is within 
the jet surface (see ’wide jet’ in Figure [5]), a strong cancel¬ 
lation takes place, like in the spherical case. However, if the 
jet is small and only one of the two regions is within the jet 
surface (see ’narrow jet’ in Figure[5l) then the suppression is 
less significant and the integrated polarization can be larger. 

From this qualitative discussion it is clear that an 
anisotropic pitch angle distribution can produce a high local 
polarization, but the polarization integrated over the emit¬ 
ting region will still be suppressed. If the visible region in¬ 
cludes only the local regions characterised by high polariza¬ 
tion, then the total polarization might be larger. However 
this requires a very narrow jet and/or fine tuned conditions 
on the location and extension of the jet. Another possibility 
to avoid such a cancellation is to consider either a emitting 
surface with no curvature or a patchy jet. In both cases, 
however, peculiar conditions between the pitch angle and 
the direction to the observer must be satisfied. Summaris¬ 
ing, from this qualitative discussion we conclude that high 
can be achieved only in hydrodynamically unphysical 
configurations (such as a narrow jet or a planar geometry of 
the emission region) and/or in contrived and unrealistically 
fine tuned configurations. Lorentz transformation (non con¬ 
sidered in Figure E} do not changes the general results of 
this qualitative discussion. 

In the following sections we present numerical estimates 
of the local and total polarization for different examples of 
anisotropic electron distributions and for different configu¬ 


rations of the magnetic field, and discuss the properties of 
the jet needed in order to maximise the circular polarization. 

7.2.1 Gaussian pitch angle distribution 

First, we consider a gaussian distribution Y{a) oc 
exp [^— j and a uniform magnetic field in the plane 

of the shock. The local polarization is derived from equa¬ 
tion [8] and can reach very high values due to the factor 
llY{ip)dY{<p)/dip. However, we need to account for the fact 
that the number of electrons with pitch angle a such that 
|q; — Qo| > 0 - 1 :, is negligible, and the flux from these regions 
is suppressed due to the lack of electrons with such pitch 
angles. Only electrons with la — ao| < ctc contribute to the 
emission. Suppression of the emission at |a — ao| > (Tq is 
assured by the multiplicative factor Y{ip) in equation (20] As 
a result, the radiation can reach the observer only from a 
small annular region for which the corresponding pitch an¬ 
gle is around oiq. This can be clearly seen in the example in 
Figure|6] (where we have chosen ao = 30° and CTq, = 5x 10^). 
As predicted by the qualitative discussion presented in the 
previous section, there is a region contributing with posi¬ 
tive values and a region contributing with negative values. 
In particular, for a > ao (a < ao) the local polarization 
is positive (negative). As can be seen in Figure [6] (upper 
panel) the maximum value is reached when the observer 
sees the edges of the distribution, i.e. for |a — qo| — (To, 
and its value is around l/Y{(p)dY{(p)/dtp ~ l/uo. Thus 
the maximal circular polarization from a point like region is 
Po^Nax ~ l/(o-a7e,min) (assuming aa'fc^.min > 1, Otherwise 
Po^max ~ 1). This explains why in the example in Figure [6] 
(where (Jc = 5 x 10“^), the maximum value of the local 
polarization (upper panel) is around ~ 20/7e,min. A larger 
polarization requires a narrower distribution, and implies a 
narrower extension of the two annular regions from where 
the radiation can reach the observer. The map for the flux 
weighted polarization is shown in the bottom panel and is 
the convolution of the map for the local circular polarization 
(upper panel) and the emissivity (middle panel). 

When integration over the emitting surface is per¬ 
formed, partial cancellation from contributions with oppo¬ 
site rotation directions takes place. We found from numerical 
estimates that the final total polarization is of the order of 
pcirc < and depends only slightly on the values 
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Figure 6. Polarization levels for the case of an anisotropic pitch 
angle distribution (a gaussian distribution peaked at ao = 30° 
and with width cr = 5 X 10“^) and uniform magnetic field. Up¬ 
per panel: maximum local circular polarization ^ 

Middle panel: weighted emissivity. Bottom panel: emissivity- 
weighted local circular polarization per unit surface. The pos¬ 
itive and negative contributions nearly cancel each other and 
pcirc ^ 0.3/7e,min- A higher polarization can be obtained by 
considering a jet whose position is selecting only (or mostly) the 
negative (or positive) region, like the ones represented by the 
green (pc*rc 2/7e,min) and orange (p=»rc bl'ye,min) circles. 
The white solid circle shows the region located at 1/P from the 
observer. 


of (to and Ua (for the particular case in Figure |6] we find 
pcirc _ o.3/7e, min). For qq = 90° there is no net polariza¬ 
tion, since the map is symmetric and perfect cancellation 
takes place. Strong cancellation can be avoided only by con¬ 
sidering a jet whose position allows the observer to see only 
(or mostly) the region contributing with positive (or nega¬ 
tive) polarization. This can be obtained, for example, with 
a very small jet with opening angle 6jst < o-ajr (see Fig. [6] 
bottom panel, orange circle). Such a tiny jet is, however, hy- 
drodynamically unphysical. A wider jet may also result in a 
polarization much larger than that obtained in the spherical 
case, as long as the edge of the jet lies on the separation 
between the positive and negative region, as the one shown 
with a green circle in the lower panel of Fig. [B] This config¬ 
uration requires, however, an extremely fine tuned position 
of the jet as compared to the observer. 

Similar results can be obtained in the case of a radial 
magnetic field. In this case, the surface of equal pitch an¬ 
gles is an annular region centred around r = 0, since tp de¬ 
pends only on r (see equation I19|l . and goes from :/? = 0° to 
= 90° going outward from r = 0 to r = 1. Again we have 
two annular concentric regions with opposite signs, one cor¬ 
responding to the rising part of the pitch angle distribution, 
where the derivative is positive, and one corresponding to 
the decreasing part, where the derivative is negative. Con¬ 
tributions from these two regions cancel themselves and the 
total polarization is still of the order of l/ye. Again, we can 
avoid cancellation only if the jet contains one region but not 
the other one, which can be realised only if extremely fine 
tuned conditions are satisfied. 

We conclude that a gaussian distribution can not pro¬ 
duce a value of in excess of l/ 7 e under realistic condi¬ 
tions. 

7.2.2 General distribution 

While in the previous section we have presented numerical 
results for the specific case of a gaussian distribution, in 
this section we consider a more general distribution, and 
present order of magnitude estimates of the total circular 
polarization. We show that the numerical results derived 
in the previous section for the specific case of a gaussian 
distribution also apply to more general functions Y{a). 

We consider a pitch angle distribution Y (a) (shown in 
figure [T] top panel) where the rising and decreasing part 
can have different widths (cri and era, respectively) and are 
separated by a region of width (T 2 where the distribution is 
flat (i.e., the number of electrons per solid angle is constant). 
An example of the polarization map corresponding to this 
distribution is shown in the bottom panel of figured for the 
case of a radial magnetic field. 

Neglecting numerical factors order unity, the local cir¬ 
cular polarization is such that P°*’'° 7 e ~ rpl/ai/a in the 
rising and decreasing part of the distribution (region 1 and 
3 respectively), and P°*’'° 7 e ~ — 1 in the flat part (region 2). 
We assume that both cri and as are much smaller than 1, so 
that in these two regions the local polarization is large. The 
total contribution coming from each region is then given 
by the local polarization multiplied by the area: ASi/ai, 
AS 2 , and ASs/as, respectively for the three different re¬ 
gions. The gaussian distribution discussed in the previous 
section can be recovered for CT 2 = 0 and cti = as. In this 
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Figure 7. The lower panel shows the circular polarization levels 
for the pitch angle distribution Y (o) depicted in the upper panel. 
The coordinate system {x,y) = {r cos (p,r sin tp) adopted to build 
the maps is shown in Fig. |2] Region tri corresponds to the inner 
annular region characterised by large and negative local polar¬ 
ization. Region 0-3 corresponds to the outer annular region with 
large and positive polarization. The two are separated by region 
(T 2 . The dashed circle corresponds to r = 1. 


case, ~ 1 /cr x {AS3 - ASi)/(AS3 + ASi). The sup¬ 

pression coming from integration over the two regions is of 
the order of cr, so that the total polarization is l/ 7 e. 

In fact, the contribution from the two regions is of the same 
order of magnitude (but with opposite signs) and almost 
cancel each other. More specifically, the sum of the two ar¬ 
eas is of the order of r^a, but their difference is of the order of 
rocr^. The final result ^ l/ 7 e is valid both for ai ~ <J 3 

and for ai <C as (or viceversa). In the last case, indeed, the 
local polarization in region 1 is much higher than the local 
polarization in region 3, but the area is much smaller, by the 
same factor. Again, the two contributions are of the same 
order of magnitude and cancel each other. 

One can wonder if it is possible to avoid cancellation be¬ 
tween the two terms by having APi << AS3, but ai > 0 - 3 , 
so that the contributions from the two regions are now very 
different. This configuration can be achieved by locating 
the two regions far from each other. This means to con¬ 
sider an intermediate region where the distribution is flat 
over a wide range of pitch angles: a 2 cri/ 3 . In this case 


ASi/ai <C ASs/as and can be neglected. The two regions 
characterised by large opposite polarization have now very 
different areas and do not cancel each other (see for example 
figure [T]). However, the total flux is dominated by region 2, 
and P^"" ~ l/ 7 e(Z\ 5 ' 3 /cr 3 - AS 2 )/AS 2 ~ l/ 7 e. 

We have demonstrated that the total integrated polar¬ 
ization is always of the order of 1 / 7 e, no matter the presence 
of regions characterised by a large local polarization. As al¬ 
ready discussed in the previous section, a possibility to avoid 
suppression is to have a jetted outflow whose extension and 
location allows to observer to see only photons with (large) 
positive or negative polarization, i.e. to see only region 1 or 
3. Let’s call cr <C 1 the width of this region. If OjetT < a 
then We stress that this result is valid only if 

Ojetr < cr <C 1 , and for a fine tuned location of the jet with 
respect to the observer. 


8 DISCUSSION: THE CIRCULAR 

POLARIZATION IN GRB 121024A 

Polarization measurements in the afterglow of 
GRB 121024A revealed that 0.15 days after the burst 
the optical radiation was circularly polarised at a level of 
0.6%, and linearly polar ised at a level of 4% , implying a 
ratio p“’' 7 P*™ ~ 0.15 dWiersema et al.ll2014ll . The claim 
of the detection was accompanied by the claim that these 
numbers are much in excess of the theoretically expected 
values of l/ 7 e, since the estimated l/ 7 e for electrons radiat¬ 
ing in the optical at the time of the detection is l/ye ~ 10 "^. 
Anisotropies in the pitch angle dis tribution have then b een 
invoked as a possible explanation (IWiersema et al.ll2014l L 

Our studies imply that, no matter the pitch angle distri¬ 
bution, if optically thin synchrotron radiation is dominating 
the emission, values of much in excess of l/ 7 e can 

be achieved only in very contrived and/or unphysical ge¬ 
ometries, that require extremely fine tuned conditions. We 
are then left with the possibility that 7 e is small, so that 
pcT-rc ^ 1 / 7 ^ can be large enough to explain the observa¬ 
tions. We then revisit the estimate of 7 e to understand if 
this is a viable possibility. 

The minimum Lorentz factor of those electrons that 
mainly contribute to the observed optical frequency i/abs = 
4.6 X 10^"* Hz is 'ye,min = y/2nmec{l + z)u°^’>/{“iPEqe). Us¬ 
ing the self-similar Blandford & McKee solution to estimate 
P and using B = •\/(327resnmpC^P^) (with n = no = const 
for a homogeneous medium and n = 3 x for a 

wind-like medium) we find respectively: 

'ye,min = 3 X 10^E"c"(ISM), (30) 

-fe,min = 1.3 X (Wind). (31) 

In order to explain the observed circular polarization, ac¬ 
cording to our numerical results, a value 'ye,min < 100 is 
needed, implying: 

E]ll2ny\^i^ > 30 (ISM), (32) 

> 13 (Wind). (33) 

Even for a large es ~ 1, the first condition can hardly be 
satisfied and the ISM case is ruled out, while in the wind¬ 
like case it is possible to account for a small 'ye,min ~ 100 
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provided that the density is large: A* > 3 x 10''. Even ad¬ 
mitting the possibility to have such a strong wind from the 
precursor star, this large density would imply a transition 
to non-relativistic velocity at early times, when the medium 
is still optically thick and no radiation can be emitted. 

In the context of optically thin synchrotron radiation, 
the possibility to explain r../ 0.6% in GRB 121024A 

as originated by electrons with small Lorentz factor is then 
ruled out. 


9 CONCLUSIONS 

In this work we have estimated the level of circular polar¬ 
ization of optically thin synchrotron radiation from a spher¬ 
ically expanding extended source moving at relativistic ve¬ 
locity. The equations for the estimate of circular polariza¬ 
tion in synchrotron radiation from a point -like region in the 
frame of the flu i d hav e been derived by iLegg fc Westfoldl 
(Il968tl . ISazono^^ lll969l l. and ISazonovI il 1972l l . However. in 
order to apply these equations to the study of emission from 
relativistic sources it is necessary i) to account for transfor¬ 
mations from the fluid frame to the observer frame, where 
the emitting fluid is moving at relativistic velocity, and ii) 
to perform integration over the whole visible emitting re¬ 
gion. We have considered these two effects and presented 
the equations that should be applied to derive from a 

relativistic unresolved source. 

We have considered three different configurations of the 
magnetic field: i) uniform in the plane of the shock, ii) ra¬ 
dial, iii) random in the plane of the shock. In this latter 
case, vanishes in each point of the emitting region. 

In a uniform magnetic field, even if the local polarization 
is of the order of l/7e, integration over the emitting region 
plays a fundamental role: for spherical outflow (or for wide 
jet or smaller jets seen on-axis) the symmetry of the sys¬ 
tem suppresses that vanishes due to the cancellation 

between regions with different directions of the polarization 
rotation. Values different from zero can be obtained only 
when the observer sees the edges of the jet and the jet is 
off-axis. In this case, values up to l/7e can be recovered. In 
a radial magnetic field, instead, the region inside the cone 
1/P contributes with negative polarization, while the region 
outside (which however gives a negligible contribution to the 
observed flux) contributes with positive polarization. In this 
case, the total value differs from zero also for the spheri¬ 
cal/wide jet case. However, also in this case its value does 
not exceed the limit l/7e. To summarize, for isotropic pitch 
angle distributions the value of p^"’^ is limited to be smaller 
than l/7e, and Lorentz transformations and integration over 
the surface might reduce the observed level of polarization, 
that is always similar or smaller than the one estimated for 
a point-like source observed in its comoving frame. 

For all the considered configurations we have also de¬ 
rived the value of the ratio p^"’^ j p*™. Contrary to what is 
generally stated, the ratio between circular and linear po¬ 
larization is not necessarily of the order of l/7e, but can 
indeed assume any value, depending on the configuration 
of the magnetic field and on the geometry of the emitting 
region. 

We have investigated the possibility to observe a high 
level of > l/ 7 e as a result of an anisotropic distri¬ 


bution of the electron pitch angles, motivated by the recent 
detection in the optical afterglow of GRB 121024A. We have 
demonstrated that invoking a large anisotropy does not help 
to overcome the limit l/7e. The only difference between 
the isotropic and anisotropic case is the following. For an 
isotropic distribution, the circular polarization of radiation 
coming from a point-like region is strongly suppressed, due 
to average over electrons with different pitch angles. For an 
anisotropic distribution the lack of electrons with some pitch 
angles can limit the cancellation and give rise to high val¬ 
ues of local polarization. However, cancellation now arises 
from the integration over the emitting region. The numeri¬ 
cal results and approximate analytic calculations presented 
in this work show that after integration, the total circular 
polarization is still of the order of l/7e. A higher p'^®’"'® can 
be reached only by considering a jet satisfying unphysical 
conditions (as a planar emitting region or a tiny jet) and/or 
very unlikely conditions, as an extremely fine tuned location 
of the jet as compared to the observer. The application of 
our study to GRB 121024A leads us to conclude that the cir¬ 
cular polarization measured in the optical afterglow of this 
GRB cannot be explained in the context of optically thin 
synchrotron radiation, no matter the configuration of the 
magnetic field and the shape of the pitch angle distribution. 
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